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INTRODUCTION 



Coppens and Sanders [1,2] developed a non-linear acoustic 
model with the dissipative term phenomonologically describing the 
viscous and thermal energy losses actually encountered at the 
walls of a rectangular rigid cavity. Several researchers [3,4] 
have examined the problems and experimental results are in 
excellent agreement with the theory except when a degeneracy 
exists in the cavity. 

The purpose of this research was to examine the 
finite-amplitude behavior within a cavity for which perturbation 
effects can be accurately measured and to compare the results with 
the predictions of the model. 
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II . BACKGROUND AND THEORY 

In 1975 Coppens and Sanders [2] formulated a perturbation 
expansion for the non-linear acoustic wave equation with a dissi- 
pative term describing the measured absorption properties and the 
measured resonance frequencies for standing waves within a real, 
fluid filled, rigid-walled cavity. This model predicts that, when 
a near-degeneracy exists, geometrical perturbations provide a 
mechanism whereby a nearly degenerate mode can affect the 
finite-amplitude behavior of the cavity. 

The non-linear wave equation of the viscous fluid [2] is, 

(c^^Q^ + 3;t/3t) p/PqC^^= 3^/3 [ (u/c^) 

where Cq 2 = (3p/3p) (adiabatic), Pq is the equilibrium 
density of the fluid, p the acoustic pressure, u the magnitude of 
the acoustic velocity, Y the ratio of heat capacities , and <iC an 
operator describing the physical processes for absorption and 
dispersion. The term on the right can be interpreted as a 
distribution of virtual sources created by the self-interaction 
of the standing waves. 

Pressure standing waves in a rectangular cavity of dimension 
, Ly, have the form 

p cos kjfX cos kyy cos k^z cos tot 

where 



= jTT/Lx 


j = 0,1,2, .. . 


= ^■JT/Ly 


1 = 0,1,2 


m TT Li 2 


rn 0,1,2, • . . 



and j+£+m^O 
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If the cavity is being driven near resonance, the contributions of 
non-resonant terms are negligible with respect to the resonant 
terms; the acoustic field has the form 



CO 



( 2 ) 



\«^iere 




M is the Mach number \ u/cq 1 » Ri=l/ and Rj^ is the relative 
amplitude of the n-th standing wave. By substituting (2) into 
(1), a set of coupled, non-linear, transcendental equations is 
obtained , 



where N = 1 for axial, 2 for tangential, and 3 for oblique 
standing waves. 



frequencies near the resonance frequency of the (0,1,0) standing 
wave, only the family members (0,n,0) will contribute 
significantly to the finite amplitude behavior. 

In this research we were interested in this behavior when a 
degeneracy exists between the (0,2,0) and (1,0,0) modes and the 
cavity is not a perfect right parallelepiped. 

If we define a perturbation parameter ^ as a dimensionless 
measure of the magnitude of any irregularities on the cavity 




(3) 




If a geometrically-per feet rectangular cavity is driven at 
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surfaces compared to the effective dimensions of the cavity, the 
effect of the perturbation on the standing wave Pj^ that would 
exist in the ideal cavity can be expressed as a small correction 
p' so that the true pressure field Pj^' is 
Pn‘ = Pn + eP' 

For the case of interest, n denotes the (0,2,0) standing wave and 
p' is the perturbation-generated (1,0,0) standing wave of nearly 
identical resonance frequency. 



I 
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III. SAMPLE CALCULATION OF THE PRESSURE DISTRIBUTION WITH A 



WEDGE PERTURBATION 

The pressure distribution of the (0,2,0) wave has the form 

P 2 = P cos (2TTy/Ly) cos(2cot+92) (1) 

where P and ©2 are constant determined by the driving 
conditions . 

For this example, let the equation for the perturbed boundary 
(Fig. O) be 

X = L^ {l-(AAjj) [ (4/L^) (v-3Ly4) ] [U(y-3Ly4) -U(y-L^) ]} 

= L^ [I4€f(y,z)] (2) 

where 



£ H A/L 

X 



f(y,z) = -(4/L^) (y-3Ly4) [U(y-3L^4) -U(y-L^) ] 



0 if V < 0 

U(y) = 

1 if y > 0 



D 



Differentiating P 2 and f(y,z) with respect to y gives 



9p2/3y = -(2 -;tP/L^) sin(27ry/L^) cos(2cot+02) 
and 

3f/3y = -(4/L^) [U(y-3L^/4)-U(y-Ly) ] 

-(4/L^) (y-3L^/4) [5 (y-3L^/4)-5 (y-L^) ] (3) 

and •S (y) is the Dirac delta function. 

Now, it can be shown by standard perturbation analysis [4] that 
the perturbation correction p' must satisfy the boundary 
condition 

OpV3^)l = L^Of/3y) OP2/3y)l 

X X 

Substituting, we obtain 

(3p'/3x)_ = A sin(2Try/Lj cos ( 2w t+0 -, ) t U (y- 3L /4)-U(y-L ) 

j-'jj y ^ y y 

+ (y-3L^/4) [5 (y-3L^/4)-5 (y-L^) ]} (4) 

where 

A = 8-itPL /L ^ 
x' y 

Equation (4) can be expanded as a Fourier series, 

00 

(3p'/3x) = A cos(2cot+0^) L [A cos (miry /L ). +b sinCmiry/L IJ (5) 

X m=0 ^ ^ 
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Only Aq is needed to determine the first order correction term 
due to the degeneracy between the (0,2,0) and (1,0,0) waves: 



= 1/27T 

The perturbation correction p' is 
2 2 

P' = (Ac^ /irL,^4aj ) Q]_oo®^^'^100 (ttx/L^) cos (2ojt+8 2 +t^qq) 
where the phase angle ^ is given by 

“ ^100 ^ ^ 

The pressure at the microphone position x=L^ can be written as 
follows (with 2o) / C q= 47t /L y ) , 

Pmic = ?2 + =P’ 

2 

= P [cos (2cot+0 2 ) + (£/2tt ) Q^QQSinT^QQCOS (20)t+9 2+T2.00^ ^ 

The total pressure amplitude at the microphone will be 
I Pmicl “ P)/(l+B cos "loo^^ ^in ^loo)^ 

where 

B = (£/2tt2) Q]_qq sin Tj_qq 

B can also be expressed in terms of Aq for a wedge anywhere on 
the wall x=L^ by 

B = (£ Aq/tt ) Q]_oo ^100 



( 7 ) 

C8), 
(9 a) 

(9b) 



/ 



where A^is calculated for any arbitrary position and wedge 
dimensions by the same method as developed above. 

[In reference [5], the author made a error in formulating Equation 
(2); he used y instead of (y-aLj^) and got an incorrect value of 
Aq*] 
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IV. DEFINITION OF SOME PARAMETERS 



Several parameters will be defined for the purpose of 
simplifying the mathematical formulation and elucidating the 
physical content of the equations. 

A. STRENGTH PARAMETERS: S 

S = MgQ^ 

where 3 = (1/2) (1 +y) and Q]_ is the quality factor of the driven 

(0,1,0) standing wave. The strength parameter characterizes the 
strength of the finite-amplitude interaction. It is interesting 
to note that S is one half the Goldberg number. Since the 
microphone sensitivity S^, (obtained with a B&K 4220 pistonphone) 
is known, 

S = M3Qi = 7.07 X 10-3 v^Qi 

where 

M = V]_/PCq2Sj^, S =1.2 for air, 

P = 1.293 kg/m3, Cq = 345 m/ s 

= RMS output voltage of the first harmonic component, 

and 

is the quality factor of the driven (0,1,0) standing 
wave . 

B. FREQUENCY PARAMETER: F^ 

Fn(f) = QnCl-(fn/f)^^ 

where f is the driving frequency, and fj^ and are the 
resonance frequency and quality factor for the (0,n,0) wave. 

If Ff^<l, the n^^ harmonic of the driving frequency lies 
within the half-power frequencies of the resonance curve for the 
(0,n,0) wave. 
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For reasonably large values of Q^/ 

Fi (f) = 2 Qi (f-f]_)/fi 

C. HARMONICITY COEFFICIENT: E(n) 

E(n) = (fn-nfi)/nfi 

E(n) characterizes how well the modes of a given family are 
tuned (harmonic). If | E(n)| <4, the corresponding harmonic will be 
strongly excited. 
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V. APPARATUS 



A. THE RECTANGULAR CAVITY 

The cavity (Fig.l), constructed from 0.75-in. aluminum, has 
interior dimensions 12.00 in. long, 2.50 in. high, and a width 
that can be varied between 5.50 in. to 7.00 in. in 0.25 in. 
increments. All joints were sealed with a thin layer of silicon 
grease . 

Figure 2 shows the pressure distribution for several of the 
lower modes of this cavity, and Table 0 presents the theoretical 
eigen frequencies calculated for a cavity 12 x 6 x 2.5 in. Note 
that the (0,2,0) and (1,0,0,) modes are predicted to be 
degenerate. Experimentally, the resonance frequency of the 
(0,2,0) standing wave was about 3 Hz higher than that of the 
(1,0,0) standing wave. The various configurations of wedge- 
perturbations are shown in Fig. 3. 

The source piston is set flush with the bottom of the cavity 
as near to the wall y=0 as practical, and halfway between the wall 
at x=0 and x=Ljj . In this position it can efficiently excite the 
(0,n,0) family of waves without appreciably exciting the (1,0,0) 
wave. To determine the proportion of the (1,0,0) mode, an 
auxilliary driver, an (ID-30), can be inserted at Position A. 

A microphone at Position B will sense the pressure of all 
standing waves. If the microphone is placed at Position A, it 
senses the (1,0,0) wave with minimum contamination from the 
(0,2,0) wave and at Position C it senses the (0,2,0) wave with 
minimum contamination from the (1,0,0) wave. 
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B. APPARATUS 



A block diagram of the apparatus is shown in Fig. 4. A GR 
11 61 -A coherent decade frequency synthesizer is used to produce a 
driving signal precise to within +0.001 Hz. This signal is 
applied to a 2120 MB power amplifier which, in turn, drives the 
shaker. The motion of the piston was continuously sensed by an 
Endevco Model 2215 Accelerometer mounted within the piston, and 
the output of the accelerometer was observed on a Model 130BR HP 
oscilloscope and measured on a HP 400D Voltmeter. A Schlumberger 
spectrum analyzer was used to measure the harmonic distortion in 
the piston motion. 

The output of the B&K l/4*in. microphone (with matching 
preamplifier B&K 2810) was fed into three devices; (1) an HP 400D 
VTVM to measure overall voltage level, (2) a Schlumberger spectrum 
analyzer to display the spectrum of the waveform, and (3) two HP 
302A wave analyzers to measure the amplitudes of the first two 
harmonics of the pressure waveform. 
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VI. DATA COLLECTION PROCEDURE 



Since the resonance frequencies of the cavity were observed 
to vary with time, the system was allowed to warm-up for at least 
one hour prior to data collection. The piston was then driven at 
the maximum amplitude to be expected during the run, and the 
harmonic content of the accelerometer output was analyzed and the 
piston was adjusted until the second harmonic of the accelerometer 
output was at least 50 dB below that of the fundamental. Figure 5 
shows typical results for the percent second harmonic in the 
acceleration with the piston driven at a rather large amplitude. 
Finite-amplitude measurements were limited to those frequencies 
for which V 2 /Vj^< 0 . 01 . 

There are three steps necessary for collecting accurate data; 
(1) pre-run infinitesimal-amplitude measurements, (2) finite- 
amplitude measurements, and (3) post-run infinitesimal-amplitude 
measurements . 

A. PRE-RUN AND POST-RUN INFINITESIMAL-AMPLITUDE MEASUREMENTS 

During these measurements, the piston was driven with the 
a„cce lerometer output less than 0.1 V. Observation of spectra of 
the pressure waveforms obtained at these low amplitudes showed 
that the amplitudes of all overtones were at least 60 dB below 
that of the fundamental. The resonant frequencies f^^ and the 
quality factors for the first few (0,n,0) waves were 
determined from 
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fn = (fu+f£)/2 
Qn = fn/(fu-f£) 

where and are the upper and lower half-power frequencies. 

The resonance frequency and Q of the (1,0,0) wave were determined 
by the same procedure but with the cavity driven at port A. 

The harmonicity coefficients were then calculated from 

^n “ (^n -nfx)/nfx, for (0,n,0) wave 

= (f i00“^2^ /^2' (1,0,0) wave 

where f^oo resonance frequency of the (1,0,0) wave. 

The time at which each f^_J and f^ were taken was recorded. 

B. THE FINITE AMPLITUDE MEASUREMENTS 

Throughout this portion of the run, the strength parameter S 
was maintained constant by adjusting the driving voltage applied 
to the piston; with the microphone in Position B, the HP 302A wave 
analyzer was set to the driving frequency and the driving voltage 
adjusted to keep the amplitude of the fundamental of the pressure 
waveform constant as the frequency was changed. The 
frequency was increased 0.5 Hz steps through +6 Hz about the 
resonance frequency of the fundamental. At each driving 
frequency, the harmonic content of the microphone output (usually 
up to the fourth harmonic) was measured with the second HP 302A 
wave analyzer set on AFC mode. The time of each measurement was 
recorded . 
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To determine the frequency parameter at the time the harmonic 
content was measured, fj^ at the instant of the measurement of 
was estimated by interpolation of the f^^ found in the pre- 
and post-run procedures. Figure 6 shows the drift in fj^ is 
approximately linear in time. The results were presented as 
Vj^/V]_ vs F]_(f) for a given strength parameter and 
perturbation. 
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VII. RESULTS 



The results are presented in graphical form in Appendix A and 
in tabular form in Appendix B. 

Figure 7 shows the excellent agreement between theory and 
experiment obtained when there are no perturbations so that only 
one family of waves is excited. Some discrepancies were always 
observed when the driving frequency was far from the fundamental 
resonance frequency (for ] >2). In this region, the piston 

had to be driven hard to keep constant, causing significant 
second harmonic to appear in the piston waveform (as shown by 
Fig. 5), thereby introducing a linear ly -gener ated second harmonic 
into the cavity. 

For all cases studied, the pressure was calculated following 
the methodology of Sect. III. 

For the boundary perturbation of Fig. 3, the results are 
plotted in Figs. 8 through 15. The thin line is the theoretical 
prediction in the absence of any perturbation and the thick line 
includes the perturbation correction. The circles are the 
experimental results. Along the edge of each figure is the 
information about the cavity configuration. 

A. THE WEDGE AT THE CORNER OF THE CAVITY 

Figures 3 through 11 show the results for a wedge at the 
corner of the cavity. The wedge used are those of Fig. 3 (a) -(d) 
respectively. 

3. THE WEDGE AT THE CENTER OF THE LONG WALL 

Figures 12 through 15 show the results for the same wedges 
[Fig. 3(e) -(h)] but now located at the center of the long wall. 
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For Fig. 12 and 13, = 0.177 with A for Fig. 13 twice that for 

Fig. 12. For Fig. 14 and 15, Aq = 0. 500 with ^ the same as for 
Fig. 12 and 13. 
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VIII. CONCLUSIONS 



A. The experimental apparatus and procedures are capable of 
providing data sufficiently precise to verify the prediction of 
the theory in the absence of geometrical perturbation if the range 
of driving frequencies is restricted so that the absolute value of 
the frequency parameter is less than 2. 

B. To verify the correction for a geometrical perturbation, 
it is useful (1) to have a sufficiently large correction to the 
unperturbed prediction and (2) that the sign of the correction 
change for an absolute value of the frequency parameter less than 
2 . 

C. The experimental results show that geometrical 
perturbation alters the finite amplitude behavior of the cavity, 
and that the nature these changes are in qualitative agreement 
with the predictions of the theory. However, quantitative 
agreement was not observed. 

D. Sources of the difficulties in obtaining good agreement 
might be (1) the inability to experimentally satisfy conditions of 
B above, and (2) the higher order terms neglected in developing 
the theory may not all be small. For example, terms of order 
higher than first order in Eq. ( 4) of Sect. II and terms of order 
higher than and all bjj, in Eq. (5) of Sect. IV. 
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Fiq. 4. Block Diagram 



r.'T 



[ 



r 



i ■ ■■[ ■■ i 






; >’ 



i I 






! t 



1 . - 



..i 

: 1 

— j. 

• 1 



! I. 

i i 



i ! 



1 r 





















I ‘ 



: 1 

I I 






Figure 5 

Typical harmonic distortion in output of the. accelerometer. 
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Typical time dependence of the resonance frequency of 
the lowest standing wave (0,1,0). 
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TABLE 0 



Mode 


Frequency (Hz) 


Mode 


Frequency (Hz) 


010 


565 


130 


2039 . 4 


020 


1131 


140 


2529.0 


030 


1697 


150 


3032 . 0 


100 


1131 


111 


2980.5 


110 


1205 


101 


2757.5 


120 


1599.5 


Oil 


2926.4 
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TABLE I (Fig. 7) 



mode 


time 


4 : 

""u 






f 

r 


Q 


E 




FP 


010 

020 

030 

100 

'^l = 


0 

3 

6 

9 

lldB 


time 


565.77 

1131.84 

1698.27 

1131.87 

f 


562.88 564.33 

1127.79 1129.81 

1693.52 1695.90 

1127.59 1129.73 

^2 ^3 ''4 


195.40 

278.69 

356.51 

264.20 

^2/^1 


0 

0 . 10 x 10 

0.17x10 

-0.74x10 

P 3 /P 1 


,-2 

,-4 

^ 4 /. 






12 


558.5 


-50.0 


-70.5 


-70.5 


0.0113 


0.0011 


0.0011 


-4.19 






13 


559.0 


-47.9 


-72.0 


-69.0 


0.0143 


0.0009 


0.0013 


-3.85 






14 


559.5 


-46.9 


-73.8 


-67.8 


0.0160 


— 


0.0014 


-3.50 






15 


560.0 


-46.0 


-72.6 


-67.5 


0.0179 


— 


0.0015 


-3.15 






16 


560.5 


-44.8 


-71.2 


-69.7 


0.0204 


0.0011 


0.0012 


-2.81 






17 


561.0 


-43.7 


-69.7 


-73.5 


0.0233 


0.0012 


— 


-2.47 






18 


561.5 


-42.4 


-67.8 


-76.5 


0.0269 


0.0014 


— 


-2.12 






19 


562.0 


-41.0 


-65.7 


-75.9 


0.0316 


0.0018 


— 


-1.78 






20 


562.5 


-39.4 


-63.2 


-82.5 


0.0380 


0.0025 


— 


-1.43 






21 


565.0 


-37.8 


-60.3 


-85.0 


0.0457 


0.0034 


— 


-1.09 






22 


563.5 


-35.9 


-56.6 


-74.9 


0.0569 


0.0052 




-0.74 






23 


564.0 


-34.0 


-52.2 


-67.5 


0.0708 


0.0007 


0.0015 


-0.40 






24 


564.5 


-32.2 


-47.6 


-60.6 


0.0876 


0.0148 


0.0033 


-0.05 






25 


565.0 


-31.1 


-45.0 


-54.5 


0.0994 


0.0200 


0.0067 


0.29 






26 


565.5 


-31.6 


-43.9 


-52.9 


0.0933 


0.0226 


0.0081 


0.64 






27 


566.0 


-33.3 


-46.7 


-56.0 


0.0772 


0.0164 


0.0056 


0.98 






28 


566.5 


-35.3 


-51.0 


-63.2 


0.0610 


-0.0100 


0.0025 


1.33 






29 


567.0 


-37.6 


-55.5 


-70.8 


0.0484 


0.0060 


0.0010 


1.67 






30 


567.5 


-39.8 


-60.0 


-73.6 


0.0365 


0.0036 


— 


2.02 






31 


568.0 


-42.1 


-64.3 


-74.7 


0.0279 


0.0022 


— 


2.36 






32 


568.5 


-44.1 


-68.5 


-74.2 


0.0221 


0.0013 


— 


2.71 






33 


569.0 


-46.2 


-69.9 


— 


0.0174 


0.0011 




3.05 






34 


569.5 


-47.2 


-68.7 


-67.0 


0.0156 


0.0013 


0.0016 


3.40 






35 


570.0 


-48.0 


-65.4 


-67.0 


0.0141 


0.0019 


0.0010 


3.14 






36 


570.5 


-48.9 


-65.1 


-67.2 


0.0127 


0.0020 


0.0016 


4.09 






37 


571.0 


-49.8 


-65.2 


-67.5 


0.0115 


0.0019 


0.0015 


4.43 






38 


571.5 


-50.8 


-65.4 


-68.2 


0.0103 


0.0019 


0.0014 


4.78 






39 


572.0 


-51.7 


-65.8 


-68.6 


0.0093 


0.0018 


0-0013 


5.12 






40 


572.5 


-52.5 


-66.3 


-67.4 


0.0084 


0.0017 


0.0015 


3.47 






41 


573.0 


-54.7 


-65.1 


-66.7 


0.0065 


0.0020 


0.0016 


5.81 



mode 


time 


f 

u 


"l 


f 

r 


Q 


E 


010 


44 


566.28 


563.37 


564.82 


194.03 


0 


020 


47 


1132.74 


1128.67 


L130.71 


277.82 


0.94x10"^ 


030 


50 


1699.77 


1694.97 


1697.37 


353.84 


0.17x10-2 


100 


53 


1132.80 


1128.48 


1130.64 


262.03 


-O.SLxlO-'^ 
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TABLE II (Fig. 3) 



mode 


time 


f 




f 


Q 


E 










u 


L 


r 










010 


0 


568.27 


565.28 


566.77 


189.81 


0 






020 


3 


1132.37 


1127.62 


1129.99 


238.25 


-0.31x10- 


-2 




030 


6 


1701.89 


1696.79 


1699.34 


333.14 


-0.57x10- 


-3 




100 


9 


1132.44 


1127.71 


1130.07 


238.76 


0.69x10- 


-4 




V^=-lldB 


time 


f 


^2 


^3 ^4 


P/P, 


^3^^1 


VL 


FP 




12 


561.3 


-44.3 


-69.9 -71.3 


0.0216 


0.0011 




-3.95 




13 


561.8 


-43.2 


-69.2 -72.2 


0.0247 


0.0012 


— 


-3.61 




14 


562.3 


-41.7 


-68.3 -74.0 


0.0293 


0.0014 


— 


-3.26 




15 


562.8 


-40.2 


-66.8 -80.0 


0.0347 


0.0016 


— 


-2.92 




16 


563.3 


-38.5 


-65.2 


0.0422 


0.0020 


— 


-2.58 




17 


563.8 


-36.4 


-62.9 


0.0537 


0.0025 


— 


-2.24 




18 


564.3 


-34.3 


-60.3 


0.0684 


0.0034 


— 


-1.90 




19 


564.8 


-32.8 


-57.2 


0.0813 


0.0049 


— 


-1.56 




20 


565.3 


-32.5 


-55.0 -72.0 


0.0841 


0.0063 


0.0009 


-1.22 




21 


565.8 


-33.5 


-53.5 -67.8 


0.0750 


0.0075 


0.0014 


-0.87 




22 


566.3 


-35.1 


-51.8 -62.5 


0.0624 


0.0092 


0.0024 


-0.53 




23 


566.8 


-36.9 


-50.9 -60.9 


0.0507 


0.1002 


0.0032 


-0.19 




24 


567.3 


-37.9 


-51.2 -60.8 


0.0452 


0.0098 


0.0032 


0.15 




25 


567.8 


-38.6 


-52.8 -64.6 


0.0419 


0.0082 


0.0021 


0.49 




26 


568.3 


-39.2 


-55.3 -69.2 


0.0391 


0.0061 


0.0072 


0.83 




27 


568.8 


-40.3 


-58.7 -75.0 


0.0343 


0.0041 




1.78 




28 


569.3 


-41.5 


-61.7 -78.5 


0.0300 


0.0029 


— 


1.52 




29 


569.8 


-42.4 


-64.3 


0.0271 


0.0022 


— 


1.86 




30 


570.3 


-44.3 


-67.3 


0.0218 


0.0015 


— 


2.20 




31 


570.8 


-45.5 


-70.8 


0.0188 


0.0010 




2.54 




32 


571.3 


-46.9 


-75.0 


0.0160 


0.0006 


— 


2.88 




33 


571.8 


-48.1 


-76.9 


0.0140 


0.0005 


— 


3.22 




34 


572.3 


-49.3 


-75.2 -69.1 


0.0122 


0.0006 


0.0012 


3.57 




35 


572.8 


-50.6 


-73.5 -68.0 


0.0105 


0.0007 


0.0014 


3.91 




36 


573.3 


-51.0 


-71.7 -66.5 


0.0100 


0.0009 


0.0017 


4.25 




37 


573.8 


-51.8 


-71.1 -66.0 


0.0092 


0.0010 


0.0018 


4.59 




38 


574.3 


-52.5 


-70.3 -65.7 


0.0084 


0.0011 


0.0019 


4.93 




39 


574.8 


053.8 


-72.3 -66.7 


0.0073 


0.009 


0.0016 


5.27 



mode 


time 


f 

u 


^2 


f 

r 


Q 


E 


010 


42 


568.82 


565.95 


567.38 


197.49 


0 


020 


45 


1133.50 


1128.70 


1131.10 


235.40 


-0.32x10"^ 


030 


48 


1703.98 


1698.69 


1701.34 


321.25 


-0.48x10"^ 


100 


51 


1133.85 


1129.10 


1131.47 


238.00 


0.33x10“^ 
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TABLE III (Fig. 9) 



mode 


time 


f 

u 


^2 


f 

r 


Q 


E 


010 


0 


567.85 


565.01 


566.43 


199.52 


0 


020 


3 


1128.61 


1124.81 


1126.71 


296.58 


-0.54x10"^ 

-2 


030 


6 


1697.95 


1693.16 


1695.39 


354.31 


-0.23x10 

-4 


100 


9 


1128.77 


1124.87 


1126.82 


288.63 


-0.96x10 ^ 


Vj_=-lldB 


time 


" ^^2 




V4 P2/P1 


P 3 /P 1 


P 4 /P 1 



12 


561.0 


-37.6 


- 66.2 


— 


0.0468 


13 


561.5 


-36.2 


-64.8 




0.0550 


14 


562.0 


-34.3 


-62.8 


— 


0.0684 


15 


562.5 


-32.2 


-60.3 


— 


0.0876 


16 


563.0 


-30.2 


-57.5 


-76.0 


0.1096 


17 


563.5 


-29.5 


-55.6 


-73.4 


0.1189 


18 


564.0 


-30.8 


-55.2 


-74.0 


0.1023 


19 


564.5 


-33.3 


-55.5 


-74.7 


0.0772 


20 


565.0 


-35.5 


-55.8 


-71.8 


0.0599 


21 


565.5 


-37.3 


-58.0 


-70.2 


0.0484 


22 


566.0 


-38.9 


-60.5 


-71.4 


0.0403 


23 


566.5 


-40.3 


-61.5 


-76.0 


0.0345 


24 


567.0 


-41.7 


-61.7 


-80.0 


0.0293 


25 


567.5 


-42.5 


-61.3 


-81.3 


0.0266 


26 


568.0 


-42.6 


-60.8 


-80.3 


0.0263 


27 


568.5 


-43.1 


-62.2 


— 


0.0250 


28 


569.0 


-43.3 


-63.7 


— 


0.0243 


29 


569.5 


-43.8 


-65.9 


— 


0.0230 


30 


570.0 


-44.5 


- 68.0 


— 


0.0213 


31 


570.5 


-45.2 


-69.6 


— 


0.0196 


32 


571.0 


-45.6 


-70.9 


— 


0.0186 


33 


571.5 


-46.2 


-71.7 


— 


0.0174 


34 


572.0 


-46.8 


-72.7 


— 


0.0162 


35 


572.5 


-47.5 


-73.3 


— 


0.0150 


36 


573.0 


-47.5 


-73.5 


— 


0.0150 


37 


573.5 


-48.2 


-77.0 




0.0139 



0.0017 

0.0021 

0.0026 

0.0034 

0.0047 

0.0059 

0.0062 

0.0060 

0.0058 

0.0045 

0.0033 

0.0030 

0.0029 

0.0031 

0.0032 

0.0028 

0.0023 

0.0018 

0.0014 

0.0012 

0.0010 

0.0009 



0.0006 

0.0008 

0.0007 

0.0007 

0.0009 

0.0011 

0.0010 

0.0006 



-3.84 

-3.49 

-3.14 

-2.79 

-2.44 

-2.09 

-1.74 

-1.39 

-1.04 

-0.67 

-0.34 

0.00 

0.36 

0.71 

1.06 

1.41 

1.76 

2.11 

2.46 

2.81 

3.16 

3.51 

3.86 

4.21 

4.56 

4.91 



mode 


time 


f 

u 




010 


40 


567.99 


565.11 


020 


43 


1128.94 


1125.12 


030 


46 


1698.43 


1693.67 


100 


49 


1129.10 


1125.21 



f 

r 


Q 


E 


566.55 


196.79 


0 

. ,„-2 


1127.03 


295.42 


-0.54x10 

. ^-2 


1696.05 


355.76 


-0.2Lxl0 

. ^-3 


1127.15 


289.76 


0 . 11 x 10 



45 



TABLE IV (Fig. 10) 



mode 


time 


f 

u 


f 


f 

r 


Q 


E 


010 


0 


566.19 


563.20 


564.70 


189.30 


0 


020 


3 


1128.86 


1124.95 


1126.91 


287.70 


-0.2xl0~^ 


030 


6 


1696.18 


1691.38 


1693.78 


353.24 


-0.18x10"^ 


100 


9 


1128.46 


1124.62 


1126.54 


293.29 


-0.32x10"" 


V^=-lldB 


time f 


^2 


V V 

3 4 


p^/Pi 


VL 


P4/P1 ^ 



12 


558.7 


-42.8 


13 


559.2 


-41.8 


14 


559.7 


-40.8 


15 


560.2 


-39.7 


16 


560.7 


-38.3 


17 


561.2 


-36.7 


18 


561.7 


-34.9 


19 


562.2 


-33.0 


20 


562.7 


-30.7 


21 


563.2 


-28,9 


22 


563.7 


-28.5 


23 


564.2 


-30.4 


24 


564.7 


-32.2 


25 


565.2 


-33.7 


26 


565.7 


-35.6 


27 


566.2 


-37.0 


28 


566.7 


-39.6 


29 


567.2 


-42.1 


30 


567.7 


-43.7 


31 


568.2 


-44.5 


32 


568.7 


-45.1 


33 


569.2 


-45.6 


34 


569.7 


-46.1 


35 


570.2 


-46.6 


36 


570.7 


-47.0 


37 


571.2 


-47.2 


38 


571.7 


-47.7 


39 


572.2 


-48.6 



-69.8 0.0257 

-68.7 0.0288 

-67.9 0.0324 

-66.9 0.0367 

-65.5 0.0434 

-63.5 0.0519 

-61.0 0.0638 

-57.8 -76.3 0.9794 

-53.7 -70.7 0.1041 

-50.2 -65.3 0.1274 

-47.5 -61.4 0.1334 

047.2 -59.0 0.1078 

-48.5 -58.5 0.0873 

-52.2 -62.7 0.0733 

-56.4 -70.07 0.0592 

-58.4 -74.7 0.0501 

-59.7 -77.2 0.0314 

-61.5 0.0279 

-63.7 0.0232 

-65.5 0.0221 

-66.9 0.0197 

-68.0 0.0186 

-69.0 0.0176 

-69.8 0.0167 

-70.4 0.0158 

-70.6 0.0155 

-71.3 0.0146 

-63.9 0.0132 



0.0011 


— 


-4.04 


0.0013 


— 


-3.71 


0.0014 


— 


-3.37 


0.0016 




-3.04 


0.0019 


— 


-2.71 


0. 0024 


— 


-2.38 


0.0032 




-2.04 


0.0046 


— 


-1.71 


0.0074 


0.0010 


-1.38 


0.0110 


0.0019 


-1.04 


0.0150 


0.0030 


-0.71 


0.0155 


0.0040 


-0.38 


0.0133 


0.0042 


-0.05 


0.0087 


0.0026 


0.29 


0.0054 


0.0010 


0.62 


0.0043 


— 


0.95 


0.0037 


— 


1.29 


0.0030 


— 


1.62 


0.0023 




1.95 


0.0019 




2.29 


0.0016 


— 


2.62 


0.0014 


— 


2.95 


0.0013 


— 


3.29 


0.0012 


— 


3.62 


0.0011 


— 


3.95 


0.0011 


— 


4.28 


0.0010 




4.62 


0.0023 


— 


4.95 



mode 


time 


f 

u 




010 


42 


566.35 


563.33 


020 


45 


1129.43 


1125.5:; 


030 


48 


1697.04 


1692.26 


100 


51 


1129.12 


1125.28 



f 

r 


Q 


E 


564.84 


186.79 


0 

^-2 


1127.48 


1287.84 


-0.19x10 

«-4 


1694.65 


353.94 


-0.78x10 
.. , «-3 


1127.20 


294.16 


-0.24x10 
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TABLE V (Fig. 11) 



mode 


time 


f 




f 

r 


Q 


E 


010 


0 


566.90 


563.86 


565.38 


185.80 


0 


020 


3 


1127.19 


1123.28 


1125.23 


288.08 


-0.49x10“^ 


030 


6 


1695.31 


1690.47 


1692. 89 


349.27 


-0.19x10”^ 


100 


9 


1126.76 


1122.90 


1124.83 


291.41 


-0.36xl0~^ 


V^=-lldB 


time 


" ^2 


"^3 


^4 


P2/P1 P3/P1 


P4/^l 



12 


559.4 


-39.5 


-66.9 






13 


559.9 


-37.9 


-65.3 







14 


560.4 


-36.4 


-63.9 






15 


560.9 


-34.5 


-61.8 


— 




16 


561.4 


-32.4 


-59.1 


-78. 


.0 


17 


561.9 


-30.3 


-55.7 


-72, 


,7 


18 


562.4 


-28.7 


-55.6 


-68. 


,2 


19 


562.9 


-28.9 


-50.8 


-65. 


.3 


20 


563.4 


-30.9 


-50.0 


-62, 


,5 


21 


563.9 


-32.3 


-49.7 


-60. 


,3 


22 


564.4 


-34.7 


-50.9 


-61. 


,4 


23 


564.9 


-36.3 


-54.3 


-67, 


,0 


24 


565.4 


-37.5 


-59.3 


-75. 


.2 


25 


565.9 


-38.4 


-65.0 






26 


566.4 


-39.3 


-69.7 






27 


566.9 


-39.7 


-72.2 






28 


567.4 


-40.6 


-74.5 


— 




29 


567.9 


-40.7 


-75.4 


— 


— 


30 


568.4 


-41.6 


-73.5 


— 


— 


31 


568.9 


-44.0 


-71.5 


— 




32 


569.4 


-47.3 


-70.9 






33 


569.9 


-49.1 


-71.1 


— 


— 


34 


570.4 


-49.9 


-71.6 


— 




35 


570.9 


-50.0 


-71.8 


— 




36 


571.4 


-50.3 


-72.3 






37 


571.9 


-50.4 


-72.4 







38 


572.4 


-50.4 


-72.4 


— 




39 


572.9 


-51.1 


-71.6 







0.0376 


0.0016 


— 


-3.99 


0.045 


0.0019 


— 


-3.66 


0.054 


0. 0023 


— 


-3.33 


0.0672 


0.0029 


— 


-3.00 


0.0856 


0.0039 


— 


-2.66 


0.109 


0.0058 


0.0008 


-2.33 


0.131 


0.0059 


0.0014 


-2.00 


0.127 


0.0103 


0.0017 


-1.67 


0.101 


0.0112 


0.0027 


-1.34 


0.081 


0.0116 


0.0034 


-1.01 


0.065 


0.0101 


0.0030 


-0.67 


0.055 


0.0068 


0.0016 


-0.34 


0.048 


0.0038 


— 


-0.01 


0.043 


0.0020 


— 


-0.32 


0.039 


0.0012 




-0.65 


0.037 


0.0009 


— 


0.98 


0.033 






1.32 


0.030 


— 


— 


1.65 


0.030 


— 


— 


1.98 


0.022 


0.0009 


— 


2.31 


0.015 


0.0010 




2.64 


0.012 


0.0010 




2.97 


0.011 


— 


— 


3.31 


0.011 


— 


— 


3.64 


0.011 


— 


— 


3.97 


0.011 


— 




4.30 


0.011 


— 


— 


4.63 


0.010 


— 


— 


4.96 



mode 


time 


f 

u 




010 


42 


566.95 


563.96 


020 


45 


1127.31 


1123.40 


030 


48 


1695.54 


1690.73 


100 


51 


1126.98 


1123.12 



f 

r 


Q 


E 


565.46 


189.24 


0 

-2 


1125.37 


288.18 


-0.49x10 

-3 


1093.14 


351.57 


-0.27x10 

-3 


1125.05 


291.39 


-0.27x10 



TABLE VI (Fig. 12) 



mode 


time 


f 

u 




f 

r 


0 


E 


010 


0 


564.94 


562.04 


563.49 


194.44 


0 


020 


3 


1137.75 


1133.55 


1135.65 


270.78 


0.77x10"^ 


030 


6 


1699.31 


1694.34 


1696.82 


341.28 


0.28x10“^ 

-4 


100 


9 


1137.85 


1133.56 


1135.71 


264.96 


0.48x10 


V^=-lldB 


time 


" ^2 


^3 


^4 


D /p D /O 

-2/^1 


P4/P1 ^ 



12 


557.5 


-53.7 


13 


558.0 


-53.3 


14 


558.5 


-52.6 


15 


559.0 


-52.0 


16 


559.5 


-51.2 


17 


560.0 


-50.4 


18 


560.5 


-49.6 


19 


561.0 


-49.0 


20 


561.5 


-48.2 


21 


562.0 


-47.7 


22 


562.5 ■ 


-47.1 


23 


563.0 


-46.3 


24 


563.5 


-45.5 


25 


564.0 


-44.4 


26 


564.5 


-42.7 


27 


565.0 


-40.5 


28 


565.5 


-38.8 


29 


566.0 


-37.0 


30 


566.5 


-35.5 


31 


567.0 


-34.1 


32 


567.5 


-32.8 


33 


568.0 


-32.3 


34 


568.5 


-33.3 


35 


569.0 


-35.4 


36 


569.5 


-37.8 


37 


570.0 


-40.2 


38 


570.5 


-42.0 


39 


571.0 


-43.6 


40 


571.5 


-46.0 



-73.9 


-76.2 


0.0073 


-74.2 


-76.3 


0.0077 


-74.6 


-77.3 


0.0083 


-75.1 


-80.2 
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TABLE VII (Fig. 13) 
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TABLE VIII (Fig. 14) 
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TABLE VIII (Fig. 15) 
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